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ABSTRACT 

NICMOS 2 observations are crucial for constraining distances to most of the existing sample of z > 1 
SNe la. Unlike the conventional calibration programs, these observations involve long exposure times 
and low count rates. Reciprocity failure is known to exist in HgCdTe devices and a correction for 
this effect has already been implemented for high and medium count-rates. However observations at 
faint count-rates rely on extrapolations. Here instead, we provide a new zeropoint calibration directly 
applicable to faint sources. This is obtained via inter-calibration of NIC2 F110W/F160W with WFC3 
in the low count-rate regime using z ~ 1 elliptical galaxies as tertiary calibrators. These objects 
have relatively simple near-IR SEDs, uniform colors, and their extended nature gives superior signal- 
to-noise at the same count rate than would stars. The use of extended objects also allows greater 
tolerances on PSF profiles. We find ST magnitude zeropoints (after the installation of the NICMOS 
cooling system, NCS) of 25.296±0.022 for FllOW and 25.803±0.023 for F160W, both in agreement 
with the calibration extrapolated from count-rates > 1,000 times larger (25.262 and 25.799). Before 
the installation of the NCS, we find 24.843±0.025 for FllOW and 25.498±0.021 for F160W, also 
in agreement with the high-count-rate calibration (24.815 and 25.470). We also check the standard 
bandpasses of WFC3 and NICMOS 2 using a range of stars and galaxies at different colors and find 
mild tension for WFC3, limiting the accuracy of the zeropoints. To avoid human bias, our cross¬ 
calibration was “blinded” in that the fitted zeropoint differences were hidden until the analysis was 
finalized. 

Subject headings: supernovae: general, techniques: photometric 
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1. INTRODUCTION 

With the installation in 1997 of the Near Infrared 
Camera and Multi-Object Spectrometer (NICMOS) in¬ 
strument, the Hubble Space Telesc ope (HST) first 
gained powerful nea r-IR capabilities (iThompsonl 119921 : 
iViana fc et al.l l2009f) . With low sky and diffraction- 
limited imaging, NICMOS was ^ 10 times faster at J 
and H point-source imaging than large ground-based 
telescopes with adaptive optics. Three cameras were 
available (NICl, NIC2, and NIC3), each 256x256 pixels, 
with pixel sizes of 0"043 (NICl), 0"075 (NIC2), and 0"2 
(NIC3, which also had grism spectroscopy). The instru¬ 
ment was originally cooled to 6 IK by a block of nitrogen 
ice until lack of coolant stopped operations in 1999. In 
2002, a servicing mission installed a cryocooler (the NIC¬ 
MOS Cooling System, NCS), allowing consumable-free 
operations at 77K. 

NICMOS enabled the first probes of the earliest half of 
the expans i on history of the universe (iRie ss et alJ I2001L 
[Mil [2?in^ ISnzuki et alllWa IRubin et alJl201.3li . Al¬ 
though precis ion ground-based z > 1 SN measurements 
are possible (iTonrv et al.ll2003t lAmanullah et al.l 120101 
[Suzuki et al.ll2012L Rubin et ah, in prep), the required 
long exposure times with lOm-class telescopes make 
building a large sample expensive. NICMOS allowed for 
the measurement of precision colors (and thus, distances) 
for these distant SNe, sampling the rest-frame B, V, R, or 
I band, depending on filter and redshift. Even with the 
forthcoming Wide Field Camera 3 (WFC 3)-observed SNe 
(jCraur et ai.ll20ll iRodnev et al.l 1201411 (from the Clus- 
ter Tensing And Supe rnova search with Hubble, CLASH: 
[Postman et al.l I2012L and the Cosmic Assembly Near- 
infrared Deep Extragalactic Legacy Survey. CANDELS: 
iCrogin et al. 120111 : iKoekemoer et al.l 1201111 . NICMOS- 
observed SNe la will continue to make up the bulk of 
the z > 1 sample. 

NICMOS h as proven to b e a c hallenging instrument 
to calibrate. iBohlin et al.l (|2005[ 1 first found evidence 
of a count-rate non-linearity (CRNL) in NIC3 when ex¬ 
tending spectrophotometric standards into the near IR. 
The Space Telescope Imaging Spectrograph (STIS) and 
NICMOS showed clear disagreement over the wavelength 
range in common (8,000 to 10,000A), with NIC3 indi¬ 
cating a relative deficit of flux for fainter sources. Pa¬ 
rameterizing the CRNL in terms of relative magnitude 
deficit per dex (factor of ten in count rate), NIC3 showed 
an increase of 0.06 mag/dex for count rates from ^ 2 
to ^ 3,000/s (^ 0.18 magnitudes over this ^ 3 dex 
range). Spectroscopy and imaging from the HST Ad¬ 
vanced Camera for Surveys (ACS) agreed with STIS, 
pointing to NIC3 as the root of the problem. A com¬ 
parison of three white dwarfs against models showed a 
strong wavelength dependence to the CRNL, with the 
C RNL consistent with zero longward of 16,000A. 

iMobasher fc Riesi ([200511 first investigated this effect 
with ground-based data using both stars and galaxies. 
The stars had been observed in both FllOW (a broad 
filter spanning Y and J centered at l.l/rm) and F160W 
(similar to H, centered at 1.6/rm) in NIC2 and with 
ground-based J and H over the J magnitude range 8-17 
([Stephens et al.l [200011 . Their galaxies ranged in bright¬ 
ness down to the sky level, and were likewise observed in 
J and H, but the NICMOS data came from NIC3 instead 


of NIC2. The star measurements showed no significant 
CRNL in either NIC2 band, with the FllOW CRNL con¬ 
strained to be a factor of at l e ast 2- 3 smaller than the 
NIC3 result from IBohlin et all ([20051 1 . The galaxy mea¬ 
surements showed no significant CRNL until the mea¬ 
surements approached the sky level (J ^ 23) when the 
scatter became large and offsets ^0.1 magnitudes may 
have been indicated. The authors suggest that charge 
trapping may be responsible for the observed CRNL: ex¬ 
posures > 155s (the persistence timescale), used to mea¬ 
sure faint objects, may be able to fully fill the traps, 
r esulting in a smaller C RNL. 

Ide Jong et al.l ([200611 used exposures of star helds with 
and without counts enhanced by a flatheld lamp to di¬ 
rectly measure the linearity of NICl and NIC2. Only 
count rates between ^ 50 and 2,000 counts per sec¬ 
ond were probed by this technique in NIC2 FllOW, 
but the CRNL again seemed to be roughly constant 
in mag/dex over this range. Interestingly, the NIC2 
FllOW CRNL seemed to be the same size before and 
after the installation of the NCS and the associated 
change in temperature. In conflict with the exposure- 
time/charge-trap hypothesis, the observed CRNL is the 
same size whether the lamp-off data are taken after the 
lamp-on data (when t he charge tra p s shou ld be full) or 
before. In a dditio n, IBohlin et al.l (1200611 checked the 
IBohlin et al.l ([200511 analysis using longer grism expo¬ 
sure times. They also find the same size NIC3 CRNL 
as with the shorte r expo sures, again at odds with the 
IMobasher fc Ries^ ([2005[ 1 results and a simple picture 
of charge-trapping. (We do however note that more- 
detailed models of charge-trapping d o seem to fit lab- 
measured data, see lRegan et al.][2012ll. _ 

T aking the measurern e nts from Ide Jong et all ([2006f l 
and IBohlin et al.l ([20061 1. Ide Jond ([200611 introduced a 
routine, rnlincor, that corrects the values in an im¬ 
age using an assumed power-law relation between the 
corrected and original values. The power-law is parame¬ 
terized in units of mag/dex in the sense that 


l/{0.063[ma 

'-^P-estimated — '-^J^observed 

l/{0.029[ma 

C^rVestimated — '-"J^observed 


g/dex]/2.5-1-1} 
g/dex]/2.5-1-1} 


for FllOW 
for F160W. 


The current convention is to then use the corrected 
count rate in combination with the zeropoint provided 
from bright standard stars. This procedure was used 
to calibrate the SNe in the G reat Observatorie s Origins 
Deep Survey (GOODS) fields ([Riess et al.ll2007l l and the 
S upernova Cosmolog y Project (SCP) high-redshift SNe 
in iNobili et al.l (1200911 . 

However, this solution was not an adequate calibration. 
The 0.006 mag/dex uncertainty on the NIC2 FllOW 
CRNL translates into a ^ 0.024 mag uncertainty over the 
4 dex range between the standard stars and high-redshift 
SNe. The effect of the strong wavelength-dependence of 
the CRNL over the FllOW filter is hard to model for 
faint sources, as the amplifier glow is not at the same 
effective wavelength as the observations, and the dark 
current has no wavelength. As these sources are a sig¬ 
nificant fraction of the total background, this introduces 
^ 0.02 magnitudes of uncertainty. It is also unclear even 
what effective count rate faint observations are taken at, 
as the amplifier glow may be constant, or produced in 
short bursts of high counts/second. We note that the 
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iMobasher fc Riesi (l2005h results could indicate that the 
NIC3 FllOW power law breaks down at low count-rates 
and is wrong by ^ 0.1 magnitudes at low count-rates 
(possibly the sum of the above effects). 

Given these issues, we were awarded 14 orbitt[3 to 
complete a precision calibration of N1C2 FllOW at 
low count-rates, unlocking the full pote n tial o f the 
high-redshift SN la data. iSuzuki et al.l 1)20121 1 and 
iRubin et al.l ()2013[1 relied on a first-round SCP FllOW 
calibration against a combination of ACS WFC and deep 
ground-based J and K data. This calibration indicated a 
zeropoint 0.055 magnitudes fainter (larger) than the ex¬ 
trapolation of the higher-count-rate calibrations, show¬ 
ing a weakening of the CRNL at low count-rates. Here, 
we derive an updated result, taking advantage of the 
similar WFC3 IR bandpasses. Given the larger num¬ 
ber of archival WFC3 and NICMOS F160W observa¬ 
tions now available, we are also able derive a result for 
F160W. Using fortuitous archival observations of mid- 
redshift galaxy clusters, we make the same measurement 
for pre-NCS observations. 

Con currently with our N1C2 investigations, iRiessI 
()2010ll compared WFC3 IR starfield data against ACS 
F850LP and NIC2 FllOW and F160W. With good pre¬ 
cision (but only for count rates that are more than 10 
times higher than high-redshift SN count rates) WFC3 
IR showed a small power-law index CRNL (~ 6 times 
smaller than for NIC2 FllOW), that is approximately 
constant in mag/dex (as a function of count-rate) when 
compared against ACS and rnlincor-corrected N1C2 
images. Similar WFC3 IR CRNL me a surem ents were 
made by iRiess fc Petrol (|2010f ): IRiessI l)2011l l indepen¬ 
dently of NICMOS, so rnlincor seems to be accurate 
within the given uncertainties at these count rates. 

Figure [1] summarizes the measurements we reference. 
The size of the CRNL is shown (left axis), plotted 
against the range of count-rates over which it was mea¬ 
sured. On the right axis, we use the Union2.1 super¬ 
nova compilation, combined with BAO, C MB, and He, 
meas urements (described in more detail in ISuzuki et al.l 
l2012fl to convert from CRNL size to cosmological im- 
pact. For evaluating this impact, we use th e WQ-Wa model 
(jChevallier fc Polarskil 120011 : iLinded [200311 in which the 
equation of state parameter of dark energy smoothly 
varies with time as w{a) = wq + ■i«a(l — «)• The high- 
redshift supernovae are particularly useful in constrain¬ 
ing the time-variation, so we judge the impact using 
shifts in the best-fit Wa- A full cosmological analysis 
will be presented with other improvements in a future 
paper; for now, we compute the linear response of Wa to 
the calibration and display that linear scale. The range 
of calibrations referenced here span ~ 1.5 in Wa- This is 
twice the size of all the other statistical and systematic 
uncertainties combined. 

As the rnlincor power-law count-rate correction 
seems to be accurate at high count-rates, our strategy 
was to begin by correcting the NICMOS data for this re¬ 
lation. As all of our data (described below) encompasses 
a relatively narrow range in count-rates (centered around 
the count-rates of high-redshift SNe), we choose to derive 
an effective set of zeropoint differences between NIC2 and 
WFC3 (four, for F110W/F160W and pre-NCS/NCS). 
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Fig. 1.— Visual summary of the referenced NICMOS FllOW cal¬ 
ibration results and their approximate cosmological implications. 
Although this work is concerned with NIC2, we include NIC3 re¬ 
sults (red dotted lines) to establish the level of uncertainty in the 
behavior of the count-rate non-linearity (although these results are 
not taken at quite the same effective wavelengths as the NIC2 re¬ 
sults). Previous NIC2 results are color-coded in blue. Each line 
indicates the measured CRNL index and the range of count-rates 
it was measured at. The “first-round” result indicated a fainter 
NIC2 FllOW zeropoint at low count-rates than the other calibra¬ 
tions, which we plot here assuming that the CRNL has a constant 
size for all count-rates. T he results of our new calibrat i on are con¬ 
sistent with the results of ide Jong et all 1120061) : [Ries iOQTg) and 
are plotted in black. The cosmological results shown on the right 
axis are evaluated by fitting a t ime-varying wp-Wg model to the 
Union2.1 supernova compilation IlSuzuki et al.|[S)12l') . and aligning 
Wa = 0 with our calibration. 

This strategy captures the relevant low-count-rate cali¬ 
bration, without necessitating the interpretation of data 
in other count-rate and exposure-time regimes. 

2. DATA 

Figure [2] presents the unnormalized HST bandpasses 
referenced in this analysis. The NIC2 and WFC3 FllOW 
bandpasses are quite similar. The NIC2 F160W band¬ 
pass extends redder than the WFC3 F160W bandpass 
and requires a mild extrapolation outside the wavelength 
range of WFC3. In both cases, there is enough overlap 
that simple color-color relations can be used to cross- 
calibrate NICMOS and WFC3. For the FllOW cal¬ 
ibration, we use the F775W—FllOW color as the ab¬ 
scissa, except when F775W is not available and we use 
F814W—FllOW. For simplicity, we avoid F850LP data, 
as CCD scattering m akes the PSF quite color-dependent 
(ISirianni et al.l [T998fl . For the F160W calibration, we 
use the WFC3 F125W—F160W color as the abscissa, 
unless F125W is not available, in which case we use 
F814W-F160W or F110W-F160W. Example galaxy 
color-color relations at a range of redshifts are shown 
in Figure [S] 

As illustrated in Figure [21 the elliptical galaxy tem¬ 
plates used in this analysis are relatively flat in f\ in¬ 
side the FllOW and F160W bandpasses. We therefore 
choose to conduct our cross-calibration using Space Tele- 
sco pe (ST) magnitudes (m agnitudes that are flat in f\, 
see lKoornneef et ^ 1198611 . Selecting a different magni¬ 
tude system (e.g., AB or Vega, both of which use bluer 
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Fig. 2.— The filter bandpasses referenced in this analysis, plot¬ 
ted against wavelength. Left to right in the top panel are the ACS 
WFC F775W filter (thin solid line), ACS WFC F814W filter (dot¬ 
ted line), the NIC2 FllOW filter (filled), and the WFC3 FllOW 
filter (thick solid line). Left to right in the bottom panel are the 
WFC3 F125W filter (thin solid line), the WFC3 F160W filter (thick 
solid line), and the NIC2 F160W filter (filled). For reference, an 
elliptical galaxy template redshifted to z = 1.2 is overplotted in 
red. All normalizations are arbitrary. 


references than ST) would have resulted in different cali¬ 
bration offsets and different correlations between calibra¬ 
tion offsets and bandpass uncertainties. However, any 
cosmological results (using those cross-calibrations and 
covariance matrices) would be the same. ST magnitudes 
have the convenient advantage that the correlations can 
essentially be neglected. 

We selected our calibration galaxies from ACS images, 
looking by eye for early-type morphologies and uniform 
colors. Each of the galaxies we selected showed sta¬ 
ble colors when using photometry with different radius 
ranges (see Section |Aj2]). Stacking the ACS data for each 
galaxy and removing an azimuthally symmetric galaxy 
model (one allowed to have ellipticity and an arbitrary 
spline radial profile) revealed spiral structure in some 
galaxies; these galaxies were removed from this analy¬ 
sis. For 14 out of 28 galaxies in the F160W calibration, 
we found archival WFC3 G141 spectroscopy (covering 
IIOOOA to 17000A), allowing us to examine the near-IR 
SED and determine the redshift. (Many redshifts also 
came from the literature, as summarized in Table [TJ) 
For FllOW, where the scatter of the color-color relation 
is smaller (and thus robust if a redshift is incorreclFl) , we 
selected re d-sequence galaxies (presented for the z > 1 
clusters in iMevers et al.l 120121) from the gala xy clus¬ 
ters ISCS J1 434.4-b3426 (IBrodwin et al.l [2^ . RDCS 
J1252.9-2927 (IRosati et all 120041) . X MI^ J2235 3-255 7 
(|Mullis et al.ll2005l) . and Abell 1835 (|Abell et al.lll989l) . 
Images of the selected galaxies are shown in Figure HI 
Among our calibration galaxies are the host galaxies 

In fact, simply assuming all FllOW galaxies in the post-NCS 
calibration are at redshift 1.2 only changes the derived calibration 
by 1 milli-magnitude (mmag, 0.001 magnitudes). 



F775W - FI low (ST Magnitudes) 



F125W - F160W (ST Magnitudes) 


Fig. 3.— Galaxy color-color relations for cross-calibrating FllOW 
(top panel) and F160W (bottom panel), shown at a range of red- 
shifts. Circular markers a r e stel lar population synthesis mod¬ 
els from f^uzual &: Chariot 1120031) . while triangular markers are 
measu rements of nearby galaxies (of all types) from fHrown et aH 
II2013I) . The vertical gray lines indicate the median colors of the 
galaxies used in our calibration (for which these colors are avail¬ 
able) . 

for the SNe SCP06C0, S CPn6H5 (ISuzuk) et a].l[Wl . 
05Lan, 04Tha, and 05Red (jR.iess et al.ll2007t . For the su¬ 
pernovae blended with their host galaxies, we used only 
the supernova-free reference images in this analysis. 
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TABLE 1 

Galaxies used in this measurement 


Galaxy 

RA 

DEC 

PIDs 

Redshift 

Redshift Source 

Emission^ 

MW E(B - V) 


FllOW, Post-NCS 

FllOW.Ol 

193.22757 

-29.45461 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W_02 

193.22703 

-29.45479 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.03 

193.22706 

-29.45644 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.04 

193.23039 

-29.45358 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.05 

193.22661 

-29.45602 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.06 

193.22575 

-29.45325 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.07 

193.22816 

-29.45401 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.08 

193.23039 

-29.45451 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.09 

193.22538 

-29.45500 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

FllOW.lO 

193.22505 

-29.45262 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.il 

193.22749 

-29.45127 

s, r, c 

1.24 

RDCS J1252.9-2927 


0.075 

F110W.12 

210.27313 

2.87484 

p, m 

0.25 

Abell 1835 


0.029 

F110W.13 

210.27098 

2.86989 

p, m 

0.25 

Abell 1835 


0.029 

F110W.14 

187.35722 

1.84900 

s, j 

1.09 

Santos et ahl (1200911: 


0.022 





Dawson et al.l 112009ll 



F110W.15 

218.62519 

34.44785 

s, j 

1.24 

ISCS J1434.4-I-3426 


0.018 

F110W.16 

218.62611 

34.44568 

s, j 

1.24 

ISCS J1434.4-I-3426 


0.018 

F110W.17 

218.62549 

34.44914 

s, j 

1.23 

Dawson et al.l 1120091') 


0.018 

F110W.18 

338.83679 

-25.96012 

s, r, j 

1.39 

XMMU J2235.3-2557 


0.021 

F110W.19 

338.83591 

-25.96062 

s, r, j 

1.39 

XMMU J2235.3-2557 


0.021 

F110W.20 

338.84198 

-25.95182 

s, r, j 

1.39 

XMMU J2235.3-2557 


0.021 

F110W.21 

338.83613 

-25.96230 

s, r, j 

1.39 

XMMU J2235.3-2557 


0.021 

F110W.22 

338.83593 

-25.96250 

s, r, j 

1.39 

XMMU J2235.3-2557 


0.021 

FllOW, Pre-NCS 

FllOW.eiK.Ol 

209.95569 

62.51310 

p, b, f 

0.32 

Fisher et al.| 1199811 


0.019 

F110W.61K_02 

209.95828 

62.51344 

p, b, f 

0.32 

Fisher et al. |l998|'l 


0.019 

F110W.61K.03 

209.95741 

62.51513 

p, b, f 

0.33 

Fisher et al.l J199811 


0.019 








F160W, Post-NCS 

F160W.01 

53.07643 

-27.84864 

i, t 

1.54 

v.lSzokolv et aL| l|2004l) 

Y 

0.007 

F160W.02 

F160W.03 

53.06273 

53.06110 

-27.72659 

-27.72709 

i, o 
i, o 

1.87 

0.98 

Balestra et ahl (121^1) 

V. ILe Fevre et al.l 1200411 

Y 

Y 

0.009 

0.009 

F160W.04 

189.23715 

62.21721 

y, h, w 

1.24 

a. (Barger et al.l 1200811 

Y 

0.013 


F160W.05 

F160W.06 

F160W.07 

F160W.08 

F160W.09 

F160W.10 

F160W.il 

F160W.12 

F160W.13 

F160W.14 

F160W.15 

F160W.16 


189.23575 

189.22982 

189.25714 

189.25511 

189.03076 

189.36618 

53.15855 

53.17661 

53.16681 

53.19196 

53.18202 

53.13239 


62.21603 

62.21776 

62.20662 

62.20382 

62.16874 

62.34293 

-27.69138 

-27.69827 

-27.73859 

-27.91250 

-27.92357 

-27.81427 


y, h, w 

y, h, w 
y, h, w 

y, h, w 
w, g 

X, d 

i, o 


d, t 
1, t 


1.225 

0.95 

1.19 

1.52 

0.64 

1.15 

0.67 

0.68 

0.52 

0.73 

0.46 

0.77 


Bargeretay 12001 


Barger et al.l i 2001 


Cohen et al.l 1120001 


Barger et al.l II200I 


Wir1di_et_aBj [2^004l ' 

LeFevreet al.l 1120 


ILe Fevre et al.l 

V, [LeFevreetay _ 

vGVanzeIIa_etj!nj [2008l) 
ILe Fevre et al.l II2004I1 



0.013 

0.013 

0.012 

0.012 

0.011 

0.013 

0.009 

0.009 

0.008 

0.008 

0.007 

0.008 


F160W.17 

7.28240 

-0.93077 

k, n 

0.23 

Abazaiian et aljQ2n09|) 

0.021 

F160W.18 

137.86492 

5.85092 

z, e 

0.78 

Kneib e^^ah^^^OT 

0.045 

F160W.19 

137.86643 

5.84706 

z, e 

0.76 

Kneib 

0.045 

F160W.20 

137.86593 

5.84595 

z, e 

0.76 

Kneib 

0.045 

F160W.21 

137.86604 

5.84474 

z, e 

0.78 

Kneib e'^alT^00n|) 

0.045 


F160W, Pre-NCS 


F160W.61K_01 

137.86492 

5.85092 

a, z 

0.78 

Kneib et al. 


2000! 

0.045 

F160W.61K.02 

137.86643 

5.84706 

a, z 

0.76 

Kneib et al. 


2OTi 

0.045 

F160W.61K.03 

137.86593 

5.84595 

a, z 

0.76 

Kneib et al. 


^oogi 

0.045 

F160W.61K_04 

137.86604 

5.84474 

a, z 

0.78 

Kneib et al. 


2OT) 

0.045 

F160W.61K_05 

209.95617 

62.51328 

p, b 

0.32 

Fisher et al. 



0.019 

F160W.61K_06 

209.95872 

62.51361 

p, b 

0.32 

Fisher et al. 



0.019 

F160W.61K.07 

209.95785 

62.51530 

p, b 

0.33 

Fisher et al. 


T^'i 

0.019 


Note. — Galaxies labeled “61K” are pre-NCS. The HST Pro¬ 
gram IDs are as follows: a = GO-7887, b = GO/DD-7941, c = 
GTO/ACS-9290, d = GO-9352, e = GO-9375, f = GTO/ACS- 
9717, g = GO-9856, h = GO-10189, i = GO-10258, j = GO-10496, 
k = GO-10886, 1 = GO-11135, m = GO-11143, n = GO-11202, o 
= GO/DD-11359, p = GO-11591, q = GO-11600, r = GO/DD- 
11799, s = GO/DD-12051, t = GO-12061, u = GO-12062, v = 
GO-12177, w = GO-12443, x = GO-12444, y = GO-12445, z = 
GO-12874. 

^ This galaxy displays IR emission lines. 

2 

This galaxy displays no optical emission lines, which, at this 
redshift, likely implies no IR emission lines. 


3. CROSS-CALIBRATION PROCEDURE 

The ideal cross-calibration procedure would be to con¬ 
strain the relative amplitude of the galaxy in each hlter 
by directly modeling each pixel in each image, marginal- 
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Fig. 4.— 3” by 3” cutouts around each galaxy. The scaling is sinh so it approaches ilogarithmic at large absolute fluxes, while 
approaching linear at small fluxes. This non-linear scaling brings out faint features, such as the Einstein ring around F160W_17. 


izing out nuisance parameters for the underlying distri¬ 
bution of galaxy light on the sky, the exact alignment of 
the images, and the relative background levels. However, 
we instead selected cross-convolution for our analysis, as 
this approach limits the impact of systematics involved in 
understanding the PSF. The resulting increase in statis¬ 
tical uncertainty due to the convolution is limited by the 
convenient fact that the galaxies are significantly broader 
than the PSFs. (PSF systematics are suppressed to some 
extent when doing supernova photometry, as these sys¬ 
tematics also affect measurements of standard stars, and 
only the differential measurement is important.) 


We first resample the data onto the same pixel scale 
and orientation using astrodrizzle (jFruchter fc et aTI 
1201011 . In short, this package resamples individual expo¬ 
sures into the same (distortion-free) frame, performs an 
initial robust image combination, rejects discrepant pix¬ 
els (in the frames of the individual exposures), and then 
resamples the good pixels from each individual exposure 
to one final combined image. The name comes from the 
process of resampling, in which flux in the individual im¬ 
age is convolved with a kernel and then “drizzled” into 
a common undistorted frame. 

Using PSFs derived from bright stars, we cross- 
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convolve the images for each filter/instrument pair to be 
compared, giving the same PSF for both images (techni¬ 
cal details in Appendix lA.lll . Once each pair of images 
has the same PSF, we centroid each galaxy, then com¬ 
pute fluxes in annuli around that centroid (IA.2I1 . We 
simultaneously fit for the true radial flux of the object 
(in the cross-convolved images), the relative sky level, 
and a scaling parameter. This scaling parameter (in 
units of magnitudes) represents the instrumental color 
of the galaxy in the pair of filters considered (instrumen¬ 
tal color in that the zeropoints have not been taken into 
account). We then convert thes e instrume ntal colors to 
STmagnitude differences (IA.3I) using the iBrown et al.l 
(|2013f) galax y templates (as a cross-c heck, we use tem¬ 
plates from feruzual fc Chario^l2003D. Fin ally, we com¬ 
pute the average zeropoint offsets in IA.41 includ ing re- 
maining uncertainty in the NIC2 bandpasses fsee lA.3.f|l 
and uncertainty in the NIC2 CRNL. 

To prevent inadvertent bias towards the expected re¬ 
sults, our analysis was “blinded” (the zeropoints were 
kept hidden) until the analysis was complete. The or¬ 
der of the unblinding was as follows. First, we checked 
the code on bright standard starj^. ensuring that the 
cross-convolution code matched the results of aperture 
photometry on the input images (cal/f It/f Ic, see be¬ 
low). This is a powerful test of the PSF models, as stars 
are much sharper than the calibration galaxies. Then, we 


unblinded the F160W results, as that band is less impor¬ 
tant for the cosmological results, and could have revealed 
gross problems with the analysis. (We made no anal¬ 
ysis changes after unblinding the F160W.) Finally, we 
unblinded the FllOW. We note that only the zeropoint 
offsets were kept hidden; the dispersions were never hid¬ 
den, and provided one avenu e of f eedback for the proper 
drizzle settings (described in lA.ll) and the annuli-annuli 
correlations (described in Section iRSl) . The dispersion 
of the bright-star observations was a particularly useful 
diagnostic. 

We evaluate the systematic uncertainties by changing 
assumptions one at a time (e.g., changing the minimum 
inner annuli radius of the photometry) and rerunning 
the analysis. To be conservative, the entire range (i.e., 
maximum — minimum) is taken to be the Icr-size of that 
systematic. We sum these differences in ouadratureF^ 
The full details of the uncertainty analysis are in Ap¬ 
pendix [B]) while the contribution from each uncertainty 
is presented in Table [2l The composition of the total 
uncertainty depends on the calibration, but statistical 
uncertainty, PSF systematics, and the calibration of the 
galaxy templates share the bulk of it. WFC3 ha s its 
own calibration uncertainties, which we evaluate in IB. 81 
These are currently comparable to the cross-calibration 
uncertainties, but may be reduced with future calibration 
programs. 


TABLE 2 

Uncertainties present in the cross-calibrations. 


Uncertainty 

FllOW 

F160W 

pre-NCS FllOW 

pre-NCS F160W 

Statistical 

10 mmag 

6 mmag 

8 mmag 

8 mmag 

Calibration ot Color-Color 

1 mmag 

2 mmag 

1 mmag 

2 mmag 

Encircled Energy Correction 

2 mmag 

2 mmag 

2 mmag 

2 mmag 

PSE Shape 

8 mmag 

7 mmag 

13 mmag 

9 mmag 

NICMOS Effective Bandpass 

30A 

17A 

30A 

17A 

Annuli Correlations 

5 mmag 

1 mmag 

10 mmag 

1 mmag 

i’emplates and Extinction 

3 mmag 

13 mmag 

6 mmag 

7 mmag 

ibtal 

14 mmag 

17 mmag 

19 mmag 

14 mmag 


4. RESULTS 


TABLE 3 

The results of our measurements. 


Fit NIC2 ST Zeropoint NIC2 Low-Count- STScI ST Zeropoint 

— WFC3 ST Zeropoint Rate ST Zeropoint 


FllOW 


WFC3 Revised Bandpass —3.138 mag 25.296 

Standard Bandpass —3.145 mag 25.272 

Pre-NCS, WFC3 Revised Bandpass —3.591 mag 24.843 

Pre-NCS, Standard Bandpass —3.592 mag 24.825 


25.262 

25.262 

24.815 

24.815 


F160W 


21 HST program IDs SM2/NIC-7049, SM2/NIC-7152, 

CAL/NIC-7607, CAL/NIC-7691, CAL/NIC-7693, CAL/NIC- 

7902, CAL/NIC-7904, SM3/NIC-8983, SM3/NIC-8986, 

ENG/NIC-9324, CAL/NIC-9325, SNAP-9485, CAL/NIC-9639, 
GO-9834, CAL/NIC-9995, CAL/NIC-9997, CAL/NIC-10381, 

CAL/NIC-10454, GO-10496, GAL/NIC-10725, GAL/NIC- 

10726, CAL/NIC-11060, CAL/NIC-11061, CAL/NIC-11319, 

SM4/WFC3-11439, SM4/WFC3-11451, GO-11557, GO/DD- 

11799, GAL/WFC3-11921, GAL/WFC3-11926, GO/DD-12051, 


GAL/WFC3-12333, CAL/WFG3-12334, GAL/WFC3-12341, 

GAL/WFC3-12698, GAL/WFC3-12699, CAL/WFG3-13088, and 
GAL/WFC3-13089. 

22 This procedure is not optimal in the presence of heteroge¬ 
neous statistical and systematic uncertainties. We test our results 
by computing the RMS scatter over all analyses for each galaxy, 
adding it in quadrature to the uncertainties for each galaxy, and re¬ 
fitting the mean offset. The shifts in mean offset are only 1 mmag, 
so the heterogeneous effects are small. 
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TABLE 3 — Continued 


Fit NIC2 ST Zeropoint NIC2 Low-Count- STScI ST Zeropoint 

— WFC3 ST Zeropoint Rate ST Zeropoint 


WFC3 Revised Bandpass —2.378 mag 

Standard Bandpass —2.376 mag 

Pre-NCS, WFC3 Revised Bandpass —2.683 mag 

Pre-NCS, Standard Bandpass —2.679 mag 


25.803 

25.789 

25.498 

25.487 


25.799 

25.799 

25.470 

25.470 


Note. — The bolded items are the recommended results. The 
difference in zeropoints is the quantity described in Section lA.4l 



Galaxy Observations, 
described in this work 



Fig. 5. — Diagram for the paths between the zeropoints dis¬ 
cussed in this work. The “Blue/CALSPEC” zeropoints are ref¬ 
erenced to CALSPEC calibration stars like solar analogs and 
Vega, which are bluer than the ST magnitude ref erence (flat i n 
f\). The WFC3 CRNL m e asure ments come from [Riessl l)2010l L 
IRiess fc Petrol II2010I) ; IRiesd II2011I) . 

Table El presents our ST magnitude zeropoint differ¬ 
ences for both the standard and revised WFC3 band- 
passes (WFC3 bandpasses are discussed in IA.3.211 . We 
recommend the analyses highlighted in boldface type; 
any revisions to the WFC3 bandpasses can be interpo¬ 
lated from the pair of numbers from each result. Like¬ 
wise, any updates to the understanding of the WFC3 
zeropoints at low count-rates can be propagated through 
the results presented here into the NIC2 zeropoints. (The 
correlations between NICMOS and WFC3 zeropoints 
specified here should be taken into account in cosmo¬ 
logical fits using the SN la Hubble diagram with both 
NIC2 and WFC3-observed SNe.) 

Here, we illustrate the application of these results to 
the NIC2 zeropoints at low count-rates, ap plicab le to any 
NIC2 cal files processed with the steps in lA.ll We note 
that our zeropoints assume l"-radius encircled energy 
corrections for FllOW of 0.935 and 0.917 for F160W 


(we use PSF photometry for the supernova data, but the 
PSFs are normalized to these values), discussed further 
in lA.ll 

Figure [ 5 ] summarizes the paths for moving between the 
zeropoints we reference. Our calibration cross-calibrates 
WFC3 and N1C2, so we start with the WFC3 zero- 
points. The observed Vega WFC3 FllOW high-count- 
rate zeropoint (with our suggested bandpass revision) is 
26.072. Accounting for the WFC3 CRNL, the zeropoint 
at low count-rates is 26.032. Converting to ST magni¬ 
tude (using our bandpass revision) gives an ST zeropoint 
of 28.434. Applying our cross-calibration gives 25.296 
for N1C2 ST. This same sequence was applied to both 
FllOW and F160W; the resulting zeropoints are shown 
in second column of column of Table El For comparison, 
we take the N1C2 Vega STScl zeropoints, and convert 
to ST zeropoints using the low-count-rate conversions in 
Table El We follow this process, rather than using the 
STScl ST zeropoints, as the N1C2 Vega-to-ST conver¬ 
sion will depend on count-rate. These results are in the 
final column of Table El Our zeropoints range between 
0.004 fainter (higher) for post-NCS F160W to 0.034 mag¬ 
nitudes fainter (higher) for post-NCS FllOW, but show 
reasonable consistency. Other low-count-rate zeropoints 
(Vega or AB) can be computed using the low-count-rate 
offsets given in Table El We remind the reader that in¬ 
terpreting the photometric measurements should be done 
using a modified bandpass, as the C RNL p referentially 
affects blue wavelengths (discussed in lA.3.1|) . 

As a modest related result, we also note that the 
galaxy-galaxy scatter in the zeropoint estimates is a few 
percent. This limits spatial variation in the N1C2 CRNL 
to ^ 10 percent, at least on ~ 1" scales. 

5. CONCLUSIONS 

This work presents a cross-calibration of the 
N1C2/WFC3 FllOW and F160W zeropoints at the the 
low count-rates applicable to high-redshift SNe la obser- 
vati ons. These measurement s are in tension with both 
the iMobasher fc Riessl (|2005ll results (at least 0.1 mag 
tension), and some earlier unpublished SCP work (0.03 
magnitudes tension). We note that this tension is not 
due to the version of calnica; we get essentially the 
same NICMOS magnitudes with the improved version 
4.4.1 as with the older 4.1.1 that the pre-2008 results 
were run with (see the discussion of the improvements in 
iDahlen et al.ll2008ll . Our results show no tension with the 
higher-count-rate zeropoint and CRNL measurements, 
with our results having smaller uncertainties at low count 
rates. A new “Union” compilation of SNe using this cal¬ 
ibration will be presented in a future paper. 
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Facilities: Hubble Space Telescope. 


APPENDIX 

A. DETAILS OF THE CROSS-CALIBRATION 
A.l. Data Processing 


We started with the NICMOS cal files, flat-fielded files that have had cosmic rays rejected using the “up -the-ramp” 
multiple readouts. We processed each cal file first with the STSDAS task pedsky (IBushouse et al.1 [20001) . to remove 
the variable quadrants seen in NICMOS data. We then ran rnlincor, to correct the images for the CRNL as measured 
at high count rates. After this process ing, amplifie r glow and other forms of spatially-variable background remained, 
so we ran the subtraction detailed in iHsiao et al.l (|201Clll . We used either the “low” or “high” background models, 
selecting the one that minimized the median absolute deviation of the imagel^ We masked the erratic middle column, 
rather than attempting to recover the flux, as this data is far less important for our extended objects than for the SN 
data with which that work was concerned. Even after pedsky and our sky subtraction, the sky level in each image 
varies spatially. We thus fit for the residual sky under each galaxy, as shown in Equation I A2I 

We did no pre-drizzle processing of the WFC3 (fit, calibrated flat-fielded exposures) data or the ACS (flc, cali¬ 
brated, flat-fielded, charge-transfer inefficiency-corrected exposures) data except using tweakreg to alig n the images. 
To p revent cosmic ray hits in the ACS images from being considered objects, we aligned L.A.Cosmic (|van DokkumI 
l2001h cleaned images. (Many of the ACS F775W visits had only one exposure per set of guide stars, so we chose to 
always align the input images (fit/flc), rather than stacking all of the data with a given set of guide stars and align¬ 
ing the stacks.) To assist with the WFC3 alignment, we replaced each bad pixel with the median of the surrounding 
values. We then transferred the alignment to the original fit or flc images. For each instrument, we selected an 
optimal reference image based on depth and overlap with other images. 

We used astrodrizzle to resample all data to a common pixel scale (0"05, the native scale of ACS) and orientation 
(arbitrarily chosen to be North-up East-left) for cross-convolution. We selected a Gaussian kernel, with pixfrac=l 
(the FWHM of the kernel in the input pixel scale). In testing, the kernel settings only had a mild impact on the 
dispersion of measured magnitudes. To prevent the loss of flux in the cores of bright stars, we weight each pixel in 
the drizzling by the exposure time of the image (this keeps the Poisson-dominated pixels from being deweighted) 0 
In our processing, we included some data quality (DQ) values that are non-zero, but still indicate a reliable flux 
measurementOddly, the post-NGS NICMOS data showed a difference in fluxes before and after astrodrizzle of 

0.7%E!I 


In order to cross-convolve the images, we must have an accurate PSF for each filter. Even if we had perfect model 
PSFs for the observed pixels, drizzling the data onto a new set of pixels will broaden the PSFs, making empirical PSFs 
a necessity. To derive NICMOS and WFC3 PSFs, we downloaded P330E data (a solar-ana l og calibration sta r with 
many observations), and derived a convolution kernel that matches Tiny Tim ( Krisd[19931 iKrist et al.l 1201 ill PSFs 
to the drizzled P330E data. Although P330E is not as red as the calibration galaxies (and will thus have a slightly 
different PSF), it is the reddest standard for which a large amount of IR data exists. We discard any images that have 
non-zero DQ flags near the PSF core except the flags in the footnote. To derive ACS PSFs, we used bright, isolated 
stars selected from the fields, as there were not enough P330E images to derive PSFs. 

These PSFs must be normalized. For this purpose, we normalize with a circular aperture of l"-radius, which is large 
enough that the variation in encircled energy (EE) with object SED is a few mmag for all filters. It is also large enough 
to ensure t hat resampling t he im age does not affect the EE values. The normalization values for ACS F775W/F814W, 
taken from iSirianni et ffi] (|2005f) . are 0.955. For NICMOS, we compute the EE values using Tiny Tim (version 7.5) 
with a range of SN, galaxy, and standard-star SEDs. We use 7x oversampling (^ O'.'Ol per pixel), which matches the 
PSFs we use for SN photometry. (The EE values change coherently by ^ 0.2% if we use lOx oversampling in stead.) 
The av erage normalization values are 0.935 for FllOW and 0.917 for F160W. For WFC3, we use the values from iHartid 
(1200^ (to best match the STScI WFC3 calibration): FllOW: 0.932, F125W: 0.927, and F160W: 0.915. 


This ord er, rnlincor then sk y -subtraction, was th e opposite 
order of what I Suzuki et al.l 1120121k IRubin et al.l II2Q13I) did. The 
resulting difference in the supernova fluxes is only about 1%, and 
we will publish an update in a forthcoming paper. Our order here 
seems to improve the agreement between NIC2 and WFC3 at the 
lowest count-rates. 

For similar reasons, we also increase the minimum cosmic-ray- 
rejection threshold to 3.0/2.0 times the derivative (instead of the 
default 1.5/0.7), used with bad-pixel rejection algorithm, miiuned. 
Before drizzling, we also scale the NICMOS image uncertainties by 
a c onstant factor for ea ch image to achieve accurate uncertainties; 
see lSuzuki et all 11201211 for details. 

For NICMOS, we allowed pixels containing flags 512 (cosmic 


ray in up-the-ramp sampling), 1024 (pixel contains source), and 
2048 (signal in zeroth read). For WFC3 IR, we allowed flags 2048 
(signal in zeroth read) and 8192 (cosmic ray detected in up-the- 
ramp sampling). 

This scale is such that the post-NCS drz images had to be 
scaled by 1.007 to match the cal images. None of the other data 
showed the same effect after accounting for pixel-area variation. 
We verified using drizzlepac pixtosky .xy2rd that the difference 
was not due to an assumed plate scale change. As we perform the 
supernova photometry on the cal images, rescaling the drz images 
is the correct procedure. 
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A.2. Fitting the Instrumental Colors, k 

We centroid each galaxy in each drizzled stack by maximizing the flux inside a 0'/15 radius aperture (it makes virtually 
no difference if 0"1 is used instead). (The signal-to-noises of these galaxies are high enough that this procedure is 
not significantly biased.) We then extract annular fluxes, f, in 1-pixel-radius steps from 1 (or 3) to 10 (or 15) pixels, 
weighting each pixel by the fraction covered by the annulus. To obtain each color, we minimize the following expression: 

rT.c-i .r-Hlogld . (Al) 

r is the residual from the model: 

r = f - [10"°'^ ^F + sa\ , (A2) 


where F is the modeled flux of the galaxy in each annulus, s is the modeled sky value, a is the area of each annulus, 
and k is the modeled ratio of instrumental count-rates (measured in magnitudes). This is the count-rate ratio (as 
observed) between two filters and/or instruments, without correcting for the object SED or the zeropoints. There are 
arbitrary scaling and offset factors, which we handle by fixing s and k to zero for one filter. Although there is only 
one sky parameter, the symmetry of the annuli implies that the fit is insensitive to linear spatial variation of the sky 
(as well as a constant offset). 

C is the covariance matrix of the f values. The diagonal terms of C are: 


^ 10-0-4feF, + sa, __ 

— , “r VgW , 

gt 


(A3) 


where g is the gain of the image (ADU/electron), t is the exposure time, and Vgky is the sky variance as determined 
empirically from object-free regions of the image. The first term is the Poisson uncertainty on the count-rates of 
the galaxy, while the second represents sky noise. As every image gets resampled by astrodrizzle, then convolved 
with another PSF, and then integrated in annuli (which share fractional pixels between neighboring annuli), there are 
large off-diagonal correlations. These correlations {pij) are also found empirically from object-free regions; we then set 

Cij = Pij yJCiiCjj. 


A.3. Fitting Zeropoint Offsets, k^"^, for Each Galaxy 

After obtaining the k values (fitting out the F values and the s values), we can fit the inter-calibrations. For 
the abscissa k values, we scale out the following zeropoints: ACS F775W: 26.41699, ACS F814W: 26.79887, WFC3 
FllOW: 28.40001, WFC3 F125W: 27.9803, and WFC3 F160W: 28.1475 (these are the STScI zeropoints, with the 
WFC3 zeropoints shifte d by 0.04 magnitudes for the WFC3 CRNL, see Section IR^ . (The impact of the uncertainties 
on colors is discussed in IB(3I Note that only differences in these zeropoints are meaningful, as they are used only to 
measure the color of each galaxy.) This gives us the abscissa ST magnitude color for each galaxy in the analysis. We 
fit linear relations to the color-color relations (see Figure [3] for typical relations), using templates with abscissa ST 
magnitudes within 0.25 of each observed galaxy. (This ±0.25 magnitude cut ensures that our results are not affected 
by the fact that the relations are not quite linear. This cut is large enough such that we always have several templates 
available to derive a local relation.) We subtract these relations from the ordinal k values, producing estimates of the 
ST magnitud e difference between NICMOS and WFC3, 

We use the iBrown et al.l (j2013[l galaxy templates for our primary analysis. These templates are constructed using 
spectra and photometry of 129 nearby galaxies, with some interpolation using (mostly) stellar population synthesis 
models from iBruznal fc Char loti (|2003f i (plus dust and PAH compone nts). Although coii structed from nearby galaxies, 
they match observed color-color relations at z 0.4 (for details, see IBrown et al.ll2013ll . lending support to their use 
at higher redshift. As a cross-check, we use iBruzual fc Charlotl (|2003ll models, but do not use this as our primary 
analysis. 

Synthesizing the color-color relations requires knowledge of the bandpasses, especially of NICMOS and WFC3 
FllOW and F160W (because of the shallow slopes of the color-color relations, the other bands are less important). 
Uncertainties in these bandpasses are described below. 

A.3.1. NICMOS Effective Bandpass 

The NICMOS CRNL depends strongly on wavelength, so the effective bandpasses of NICMOS will depend on 
count-rate, as illustrated in Figure [H We measure excellent agreement between synthesized (using the 2014 March 
CALSPECE]) and measured magnitudes among G191-B2B, GD153, GD71, GRW±70 5824, WD1657±343, P041C, 
P177D, P330E, SNAP-2, VB8 (the data here are saturated in F160W), 2M0036±18, and 2M0559-14 (FllOW data 
only) using the synphot NIC2 bandpasses at high count-rates. This check limits any significant deviation from the 
standard bandpass to only the effects of the CRNL. There are no blue NIC2 medium or narrow-band filters, and no 
NIC2 grism, so we we cannot measure the change in NIC2 CRNL with wavelength. However, the CRNL (in mag/dex) 
is roughly linear with wavelength for NIC3, where it was measured in small wavelength bins using the grisms. We 
thus parameterize the effect on the NIC2 bandpasses using a function that is linear (in magnitudes) with respect to 
wavelength (i.e., an exp (A) bandpass warping function). This function is constrained to be 0.063 magnitudes/dex at 

http://www.stsci.edu/hst/observatory/crds/calspec.html 
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11 ,000A, and 0.029 magnitudes/dex at 16,000A, matching the high-count-rate Ide Jong et al.l (|2006D measurements of 
the NIC2 CRNL in the FllOW and F160W data, respectively. As we do not know the effective wavelength of amplifier 
glow (or how to treat dark current), we do not assume that this function should be evaluated with four dex (for the 
four dex separating the supernovae and standards). We instead parameterize the deviation from the high-count-rate 
bandpass in terms of the nuisance parameter /3 (see Section lA.41) . which warps the bandpasses by at 11,000A 

and 10“^5 0 029 at 16,000A. Our standard analysis conservatively assumes a Gaussian prior of 2 ± 2 on /?, so that both 
four and zero are easily accommodated. 

The sensitivity of NICMOS improved preferentially in the blue with the installation of the NCS. For the pre-NCS 
data, the bandpass must therefore be adjusted. Turning again to the NIC3 grism data (in G096L and G141L), we see 
that the pre/post-NCS sensitivity change is roughly linear with wavelength. As with the wavelength dependence of 
the CRNL, we fix the NIC2 slope with wavelength using the pre/post-NCS zeropoint change in the FllOW and the 
F160W (0.45 and 0.33 magnitude^). This lets us handle pre-NCS data with the same bandpass model, just with the 
above prior on /3 changed to 5.6 ± 2. 



Wavelength (A) 


Fig. 6.— Plots of our model of the NIC2 FllOW bandpass under different conditions, with arbitrary normalization. The blue line 
shows the bandpass taken from synphot. Observations of bright standard stars with a range of colors show consistency with th is bandpass 
at those high count-rates. In red, we show the assumed bandpass two dex fainter (/3 = 2). As described in Section lA.3.11 the CRNL 
preferentially acts at blue wavelengths, shifting the bandpass effective wavelength to the red for lower count-rates. The pre-NCS NIC2 had 
worse sensitivity at blue wavelengths, giving a further effective-wavelength shift to the red shown in green (/3 = 5.6). 


A.3.2. WFC3 Effective Bandpass 

Unlike NICMOS, the WFC3 CRNL is roughly independent of wavelength. Thus, establishing the WFC3 bandpasses 
at high count-rates is sufficient for all count-rates. (Although the galaxies in this analysis are close to zero ST color 
on average (flat in f\), knowledge of the WFC3 bandpasses is necessary to compute the ST magnitude zeropoint from 
the bluer standard stars.) As with NIC2, we check the observed and synthesized magnitudes of the standard stars 
G191-B2B, GD153, GD71, GRW-b70 5824, WDI657-b343, P04IG, P177D, P330E, SNAP-2, and KF06T2 (for F160W, 
there is also data for VB8). These stars span a smaller range of colors than the stars observed with NIG2, but strongly 
indicate that shifts of the bandpasses to the red are necessary. Goincidentally, the effective-wavelength shifts needed 
for both filters are 60A. We implement these shifts using the same smooth warping function used in Section lA.3.11 As 
we only need the bandpass for converting between the Calspec-derived zeropoints and the ST magnitude zeropoints, 
the choice of functional form for the effective-wavelength shift will only have a small effect. 

A.3.3. Synthesized High-Count-Rate Zeropoints 

In TablelH we present our Vega zeropoints derived from standard stars using I"-radius aperture photometry. (Vega is 
close in color to the average standard used in this determination; these zeropoints ca n be transformed usi ng Table 0) 

The WFC3 bright zeropoints are fainter than the STScI zeropointJ^. as noted by iNordin et al.l (|2ni4f l (who used 
PSF photometry). The WFG3 Vega zeropoints are almost independent of bandpass used, as the average color of the 
standard stars is not very dissimilar from Vega. Varying the photometry radius used can vary the zeropoints by several 
mmag, so these zeropoints are only tied to CALSPEC at the level of ^ 0.01 mag. The CALSPEC system itself also 
has uncertainty, so all of these zeropoints are most accurately defined with respect to other calibrations to that system. 

The NIGMOS zeropoints show good agreement with their STScI counterparts, with some scatter. We note that the 
NICMOS bright zeropoints are only presented for comparison to the faint zeropoints, and do not enter our analysis 
(except to constrain the pre-NCS NIC2 bandpasses). 

http://www. stsci. edu/hst/nicmos/perf ormance/photometrY/postncs Jiej'mOrds. html /, i ^ ^ ^ / , 

I- I- j I- ana http:^sci.edu/hst/nicmos/performance/photometry/pr6ncs_k6yi, 

29 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
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TABLE 4 

High-Count-Rate Vega Zeropoints 


Bandpass 

Observed Zeropoint 

STScI Zeropoint 

WFC3 FllOW, Synphot 

26.074 

26.063 

WFC3 FllOW, Suggested Revision 

26.072 


WFC3 F160W, Synphot 

24.708 

24.695 

WFC3 F160W, Suggested Revision 

24.708 


NICMOS FllOW, Synphot 

22.973 

22.964 

NICMOS FllOW, Pre-NCS 

22.500 

22.500 

NICMOS F160W, Synphot 

22.144 

22.153 

NICMOS F160W, Pre-NCS 

21.816 

21.816 


Note. — These zeropoints are computed using l”-radius aperture photometry of bright standards with the March 2014 CALSPEC 
spectra. The proposed revisions of the WEC3 bandpasses have little effect on the Vega ze ropoints, as t h e ave rage color of the standards 
is not far from Vega. We find fainter (larger) zeropoints for WFC3, in accordance with iNordin et aTI 1120141 1. Our high-flux NICMOS 
zeropoints are presented for comparison only, and do not enter our analysis. 


TABLE 5 

Synthesized Zeropoint Differences 


Bandpass 

Effective 

Wavelength 

ST - Vega 
(Mag) 

AB — Vega 
(Mag) 

WFC3 FllOW, Synphot 

11797 

2.3826 

0.7647 

WFC3 FllOW, Suggested Revision 

11857 

2.4024 

0.7728 

WFC3 F160W, Synphot 

15436 

3.4978 

1.2566 

WFC3 F160W, Suggested Revision 

15496 

3.5131 

1.2634 

NICMOS FllOW, Synphot 

11575 

2.2936 

0.7328 

NICMOS FllOW, Suggested Low-CR 

11605 

2.3036 

0.7366 

NICMOS FllOW, Pre-NCS and Low-CR 

11659 

2.3215 

0.7434 

NICMOS F160W, Synphot 

16159 

3.6474 

1.3147 

NICMOS F160W, Suggested Low-CR 

16175 

3.6515 

1.3165 

NICMOS F160W, Pre-NCS and Low-CR 

16206 

3.6588 

1.3196 


Note. — This table is intended to aid conversions among the different magnitude systems. We present the effective wavelength of 
each filter, computed for a source flat in f\. We also present ST — Vega and AB — Vega magnitude conversions. Each WEC3 result is 
presented with and without our proposed bandpass shift. We also present results using the NIC2 bandpasses at high-count rates, with 
the effects of the CRNL taken into account for low count-rates, and pre-NCS at low-count rates. 


A.4. Fitting the Global Zeropoint Dijferenees, 

Tests involving fitting a scale between images of the same galaxies in WFC3 data (with a range of spatial offsets and 
rotations) reveals a ^ 0.03 magnitude scatter. We take this as due to different pixel sampling in the undersampled 
images. The existence of this irreducible scatter implies that the statistical uncertainty is best judged (in part) using 
the observed dispersion of the scale factors about the mean. We must take into account residual uncorrected CRNL 
for both NIC2 and WFC3, as well as the partially known effective bandpass at these low count-rates. As we have 
enough data points to reliably estimate both calibration parameters and uncertainties using maximum likelihood, we 
minimize the following expression for each calibration: 


E 


[kF - ikF 




^mean SNJ 


-PCI 


Tnlincor\l2 


T? + Cr? , 


+ E + Tint) ■ 

i 


(A4) 


kF is the ST magnitude NIC2—WFC3 difference measurement for each galaxy (with measurement uncertainty erf). 
Mfnimcor jg amount of non-linearity correction rnlincor applies to each galaxy. It is thus a surrogate count-rate 
measurement, with lower count-rates giving higher corrections. is the mean rnlincor correction for the 

high-redshift SNe la near maximum, equal to 0.23 magnitudes in FllOW and 0.10 magnitudes for F160W. (7™*'““ is 
the change in kF with respect to a change in count-rates by one dex due to the estimated wavelength-dependence of 
the CRNL (with respect to ST magnitude); it is thus a measure of the color of each galaxy. (Emission lines also play 
a role, but most of the variation is due to color.) We present our measurements of these parameters in Table [6l 
The fit parameters are as follows. A:®"'’ is the zeropoint offset defined for zero ST color and a parameterizes 

any residual CRNL in either NIC2 or WFC3 (for simplicity, we assume that the WFC3 CRNL is proportional to the 
NIC2 CRNL). As described in Section [A. 3. 11 /3 is used to measure any deviation from the standard high-count-rate 
bandpasses. Finally, is a fit parameter representing irreducible variance (assumed to be the same for all galaxies in 
one band). The sum is usually over each galaxy. As there are not enough objects in the small NICMOS FoV to align 
separate NICMO S datasets, the sum ranges over these datasets, if more than one is present for a galaxy. As discussed 
in Section lA.3.1[ we take a prior of 2 ± 2 on /3 for the post-NCS NICMOS data, and 5.6 ± 2 for the pre-NCS data. For 
the post-NCS data, we do not take any prior on a, as we are testing for deviation from the predicted low-count-rate 
behavior. For the pre-NCS data (which uses many fewer objects), we assume that the WFC3 CRNL is 0.01 mag/dex. 
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and the NICMOS CRNL is adequately corrected over this narrow range of count-rates (as it seems to have been in this 
count-rate range for the post-NCS data), a is thus fixed to 0.01 mag/dex / 0.063 mag/dex = 0.1587 for the pre-NCS 
FllOW data and 0 .01/0.029 = 0.3448 for the pre-NCS F160W data (recall that the 0.063 and 0.029 come from the 
Ide long et al.l[20r)^ measurements of the NIC2 CRNL at high count-rates). 

Illustrations of the fits are shown in Figure[71 We note that for the FllOW data, it appears that the NIC2—WFC3 
zeropoint gap narrows at very low count-rates (visible as higher points towards the right in the left panels). It may 
be that rnlincor over-corrects NIC2 FllOW at these count-rates. Additional systematic uncertainty is likely called 
for when using rnlincor corrections greater than 0.25 magnitudes for NIC2 FllOW. 

There are t hree faint stars in the FllOW data, allowing us to use them as a cross-check. For these, we use the 
iPickle^ (jl998l l stellar library for the color-color relation. As with the pre-NCS data, we fix a, as we do not have 
enough objects over a large enough range of count-rates to reliably fit it. Large-aperture photometry on faint stars 
does not give high signal-to-noise, but we do find consistency with the galaxy results: fco"*" = —3.16 ± 0.04 using the 
revised WFC3 bandpass. 

As another cross-check, we fix a for the post-NCS data to investigate how fitting out uncorrected CRNL affects our 
results. The calibrations, using the modified WFC3 bandpasses, are only 4 mmag and 3 mmag (FllOW and F160W, 
respectively) different. These tests indicate that our mean galaxy count-rate is close to the mean supernova count-rate. 
As a similar cross-check, we have objects spanning enough of a color range in FllOW to unfix /? (although we now 
fix a for maximum statistical power). This results in a zeropoint difference of 10 mmag. Encouragingly, we find a /3 
measurement of 8.7 ± 6.2, more consistent than not with the need to modify the bandpass at lower count rates. 
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Fig. 7. — The top panels present the galaxy measurements for FllOW; the bottom panels show F160W. The t/-axis is always kf^, the ST 
zeropoint difference observed between NIC2 and WFC3, computed using f^rown et 112013) templates. The left panels show kf^ plotted 
against the size (in magnitudes) of the rnlincor NIC2 correction for the galaxy; this is a measure of relative galaxy surface brightness, 
with higher surface brightness galaxies towards the left. The right panels show kf'^ plotted against C'rnlmcor^ effect of the wavelength 
dependence of the CRNL on kf'^ as the count-rate changes by one dex. This is a measure of relative galaxy color, with redder galaxies 
to the right. The blue boxes (one for each galaxy) represent the range in results for different analyses (e.g., varying the outer photometry 
radius) for each point. The black points and error bars represent the mean for that galaxy, with the mean error bar including o-jnt- Each 
gray line is the fit for each variant (Section [BI . The green lines present the STScI NIC2 and WFC3 calibrations, with the WFC3 IR 
zeropoints moved 0.04 magnitudes brighter (smaller) to represent the uncorrected WFC3 CRNL. In the left plots, the green lines are not 
horizontal as the WFC3 CRNL has not been corrected, and thus the expected NIC2/WFC3 zeropoint difference changes with count rate. 
In the right plots, the green lines are not horizontal as they sh ow the (/3 = 2) NIC2 bandpass shift to the red representing the preferential 
loss of blue sensitivity due to the NIC2 CRNL fSection lA.3.1ll . 
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TABLE 6 

Derived Galaxy Quantities 


Galaxy NIGMOS CRNL Abscissa Color Abscissa Value crnlincor 


FllOW, Post-NCS 


FllOW.Ol 

0.232 

F775W 

— 

FllOW 

1.273 

-3.205 

± 

0.060 

-3.158 

± 

0.061 

0.00188 

F110W_02 

0.225 

F775W 

— 

FllOW 

1.303 

-3.196 

± 

0.031 

-3.148 

± 

0.032 

0.00194 

F110W.03 

0.225 

F775W 

— 

FllOW 

1.355 

-3.204 

± 

0.047 

-3.154 

± 

0.050 

0.00204 

F110W.04 

0.230 

F775W 

— 

FllOW 

1.144 

-3.082 

± 

0.068 

-3.040 

± 

0.069 

0.00161 

F110W.05 

0.233 

F775W 

— 

FllOW 

1.171 

-3.173 

± 

0.079 

-3.130 

± 

0.089 

0.00167 

F110W.06 

0.241 

F775W 

— 

FllOW 

1.160 

-3.160 

± 

0.098 

-3.118 

± 

0.099 

0.00165 

F110W.07 

0.267 

F775W 

— 

FllOW 

1.016 

-3.031 

± 

0.143 

-2.995 

± 

0.147 

0.00135 

F110W.08 

0.234 

F775W 

— 

FllOW 

1.448 

-3.140 

± 

0.044 

-3.087 

± 

0.044 

0.00220 

F110W.09 

0.220 

F775W 

— 

FllOW 

1.369 

-3.188 

± 

0.026 

-3.138 

± 

0.031 

0.00206 

FllOW.lO 

0.237 

F775W 

— 

FllOW 

1.140 

-3.203 

± 

0.051 

-3.161 

± 

0.053 

0.00161 

F110W.il 

0.253 

F775W 

— 

FllOW 

1.145 

-3.179 

± 

0.087 

-3.137 

± 

0.087 

0.00161 

F110W.12 

0.180 

F814W 

— 

FllOW 

-0.098 

-3.138 

± 

0.005 

-3.149 

± 

0.007 

-0.00060 

F110W.13 

0.211 

F814W 

— 

FllOW 

-0.125 

-3.141 

± 

0.016 

-3.155 

± 

0.016 

-0.00070 

F110W.14 

0.215 

F775W 

— 

FllOW 

1.070 

-3.218 

± 

0.018 

-3.187 

± 

0.019 

0.00123 

F110W.15 

0.218 

F775W 

— 

FllOW 

0.766 

-3.139 

± 

0.031 

-3.113 

± 

0.034 

0.00076 

Fiiow.ie 

0.229 

F775W 

— 

FllOW 

0.966 

-3.242 

± 

0.023 

-3.207 

± 

0.025 

0.00125 

F110W.17 

0.213 

F775W 

— 

FllOW 

1.240 

-3.157 

± 

0.018 

-3.112 

± 

0.021 

0.00178 

F110W.18 

0.242 

F775W 

— 

FllOW 

1.464 

-3.169 

± 

0.061 

-3.112 

± 

0.062 

0.00274 

FllOW.ig 

0.233 

F775W 

— 

FllOW 

1.410 

-3.167 

± 

0.031 

-3.112 

± 

0.037 

0.00263 

F110W.20 

0.253 

F775W 

— 

FllOW 

1.235 

-3.117 

± 

0.079 

-3.068 

± 

0.081 

0.00225 

F110W.21 

0.264 

F775W 

— 

FllOW 

1.172 

-3.160 

± 

0.097 

-3.113 

± 

0.098 

0.00210 

F110W.22 

0.226 

F775W 

— 

FllOW 

1.379 

-3.152 

± 

0.040 

-3.098 

± 

0.043 

0.00257 


FllOW, Pre-NCS 


F110W_61K_01 

0.169 

F775W 

- FllOW 

-0.065 

-3.599 ±0.006 

-3.609 ±0.010 

-0.00045 

F110W.61K.02 

0.193 

F775W 

- FllOW 

-0.125 

-3.567 ± 0.013 

-3.582 ±0.016 

-0.00067 

F110W.61K_03 

0.191 

F775W 

- FllOW 

-0.076 

-3.599 ±0.013 

-3.610 ±0.019 

-0.00050 


F160W, Post-NCS 


F160W.01 

0.091 

F125W 

- F160W 

-0.054 

-2.400 

± 

0.019 

-2.369 

± 

0.020 

-0.00013 

F160W.02 

0.094 

F125W 

- F160W 

-0.009 

-2.407 

± 

0.016 

-2.371 

± 

0.016 

-0.00061 

F160W.03 

0.105 

F125W 

- F160W 

-0.309 

-2.435 

± 

0.041 

-2.388 

± 

0.042 

-0.00059 

F160W.04 

0.096 

F125W 

- F160W 

-0.067 

-2.422 

± 

0.016 

-2.376 

± 

0.018 

-0.00065 

F160W.05 

0.105 

F125W 

- F160W 

-0.068 

-2.335 

± 

0.045 

-2.287 

± 

0.046 

-0.00069 

F160W.06 

0.094 

F125W 

- F160W 

-0.079 

-2.414 

± 

0.017 

-2.385 

± 

0.019 

-0.00032 

F160W.07 

0.104 

F125W 

- F160W 

-0.153 

-2.533 

± 

0.091 

-2.454 

± 

0.093 

-0.00121 

F160W.08 

0.102 

F125W 

- F160W 

-0.054 

-2.413 

± 

0.029 

-2.377 

± 

0.029 

-0.00019 

F160W.09 

0.094 

F125W 

- F160W 

-0.122 

-2.410 

± 

0.022 

-2.378 

± 

0.022 

-0.00037 

F160W.10 

0.090 

F125W 

- F160W 

-0.098 

-2.423 

± 

0.009 

-2.369 

± 

0.011 

-0.00077 

F160W.il 

0.086 

F125W 

- F160W 

-0.083 

-2.432 

± 

0.009 

-2.417 

± 

0.011 

-0.00019 

F160W.12 

0.093 

F125W 

- F160W 

-0.137 

-2.430 

± 

0.016 

-2.408 

± 

0.016 

-0.00029 

F160W.13 

0.087 

F125W 

- F160W 

-0.198 

-2.462 

± 

0.010 

-2.395 

± 

0.010 

-0.00083 

F160W.14 

0.090 

F125W 

- F160W 

-0.101 

-2.379 

± 

0.013 

-2.366 

± 

0.013 

-0.00020 

F160W.15 

0.093 

F125W 

- F160W 

-0.308 

-2.362 

± 

0.038 

-2.284 

± 

0.040 

-0.00108 

F160W.16 

0.093 

F125W 

- F160W 

-0.151 

-2.413 

± 

0.021 

-2.392 

± 

0.022 

-0.00032 

F160W.17 

0.062 

F814W 

- F160W 

-0.292 

-2.436 

± 

0.005 

-2.372 

± 

0.007 

-0.00081 

F160W.18 

0.084 

F125W 

- F160W 

-0.082 

-2.399 

± 

0.028 

-2.388 

± 

0.029 

-0.00020 

F160W.19 

0.080 

F125W 

- F160W 

-0.063 

-2.402 

± 

0.019 

-2.394 

± 

0.020 

-0.00015 

F160W.20 

0.085 

F125W 

- F160W 

-0.079 

-2.365 

± 

0.038 

-2.355 

± 

0.038 

-0.00018 

F160W.21 

0.088 

F125W 

- F160W 

-0.128 

-2.397 

± 

0.043 

-2.379 

± 

0.046 

-0.00028 


F160W, Pre-NCS 


F160W.61K_01 

0.082 

F125W 

- F160W 

-0.082 

-2.677 ± 0.044 

-2.666 ± 0.045 

-0.00020 

F160W.61K.02 

0.078 

F125W 

- F160W 

-0.063 

-2.691 ± 0.029 

-2.683 ±0.030 

-0.00015 

F160W.61K.03 

0.084 

F125W 

- F160W 

-0.079 

-2.666 ± 0.060 

-2.656 ±0.061 

-0.00018 

F160W.61K_04 

0.088 

F125W 

- F160W 

-0.128 

-2.758 ±0.074 

-2.740 ±0.078 

-0.00028 

F160W.61K_05 

0.066 

FllOW 

- F160W 

-0.167 

-2.758 ±0.004 

-2.710 ±0.006 

-0.00065 

F160W.61K_06 

0.078 

FllOW 

- F160W 

-0.197 

-2.727 ± 0.008 

-2.677 ±0.009 

-0.00069 

F160W.61K.07 

0.077 

FllOW 

- F160W 

-0.172 

-2.752 ± 0.007 

-2.704 ±0.011 

-0.00067 


Instrumental magnitude difference between NICMOS and WFC3. The uncertainty is the mean statistical uncertainty on these measurements. 
^ ST magnitude offset, computed using [Brown et al.| 1 120131) galaxy templates only. This uncertainty also includes variation due to photometry 
parameters. 
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B. DETAILS OF THE UNCERTAINTY ANALYSIS 
B.l. Statistical Uncertainty 

The statistical uncertainties in the fits of fEauation lA4D are 10 mmag in FllOW and 6 mmag in F160W (both 
post-NCS). As the likelihood is approximately Gaussian, these are computed using the Jacobian matrices with the 
covariance matrix of observations. 


B.2. PSF Uncertainty 

Our PSFs, derived from P330E, are not identical to the PSFs of the galaxies. This will lead to systematic mismatches 
between the photometry for different filters. We verify our PSFs by varying the inner radius used (either 1 or 3 pixels / 
0"05 or 0'.'15). We also vary the outer radius used (10 or 15 pixels/ O'.'S or 0"75). The range spanned by these changes 
is 8 mmag in FllOW and 2 mmag in F160W, which we take as a systematic uncertainty. We also try a fully empirical 
PSF (not relying on Tiny Tim as a first approximation). This makes a difference of only 1 mmag. 

B.3. Impact of Other Zeropoints on the Color-Color Relations 

The slopes of the color-color relations used to calibrate FllOW are ^ 0.03 mag/mag (with modest variation for 
different redshifts, templates, and abscissa colors), see Figure [3] This implies that the ~ 0.03 mag uncertainties on 
the ACS/WFC3 relative calibration (including zeropoints, encircled-energy correction, and the WFC3 CRNL) will 
contribute 1 mmag to the uncertainty on the FllOW calibration. For F160W, the slopes are ^ 0.15 mag/mag (again 
with modest variation), but the relative calibration uncertainties are smaller, as the abscissa colors generally both 
come from WFC3. We take a 2 mmag uncertainty for this relation. 

B.4. Encircled Energy Correction 

Our measurement is sensitive to the differential in encircled energy between NIC2 and WFC3. Future updates to 
the encircled energy corrections can be propagated into our results; for the moment, we take a 2 mmag uncertainty. 

B.5. Annuli Correlations 

The C matrices that we empirically determine ('Section lA.21) have large off-diagonal correlations. These correlations 
are determined using object-free regions, and thus lack (smaller-scale) variations such as those caused by focus changes 
or sub-pixel position variations of sharp cores. As an approximate way to investigate the sensitivity to the ratio of 
small-scale to large-scale correlations in the C matrices, we tried uniformly rescaling all of the off-diagonal elements 
by a range of values. These rescalings lower the dispersion in k by more than a factor of two for stellar observations 
(these point-source observations show the largest response). There is little variation in the fitted zeropoint values or 
their error bars for a broad range of scale values, from 0.97 to 0 (where 0 results in an uncorrelated matrix). These 
rescalings have an effect of a few mmags on the post-NCS results (summarized in Table [2]), which we take as systematic 
uncertainty. 


B.6. Uncertainty in Calaxy SEDs 

For the F160W bandpasses, the RMS residual from the color-color relation is 10 mmag, half of which we take as 
systematic uncertainty (to account for the fact that the average of our galaxies may not be the same as the average 
of the templates). Due to the similarities between the NIC2/WFC3 FllOW bandpasses, the scatter in the color-color 
calibration for F llOW is smaller (5 mmag ), as shown in Figure [3] We again take half this as systematic uncertainty. 
Switching to the lBruzual fc Char loti (|2003l) templates changes the zeropoints by < 1 mmag and 12 mmag (FllOW and 
F160W post-NCS, respec tively). It is a l so pos sible that our galaxies have more (or less) dust than the nearby galaxies 
used in constructing the iBrown et al.l (120131) templates. Adding 0.1 magnitude of CCM reddening (jCardelli et al.l 
Il989f) to the templates (with Ry = 3.1, so Ay = 0.31) changes the zeropoints by < 1 mmag and 2 mmag (FllOW and 
F160W post-NCS, respectively). There is also uncertainty on the Milky Way foreground extinction for each galaxy, 
but these uncertainties affect our results at a trivial level. 

B.7. AGN Variability 

It is possible that some of our calibration galaxies have AGN, allowing them to change brightness in the time span 
between the NICMOS, ACS, and WFC3 observations. However, any large variability would flag the gala xy as unstable 
as we change the inner aperture size. Small variability is possible, but would increase a-mt in Equation IA4I and so is 
already included in the statistical error bar. 


B.8. WFC3 Uncertainties 

Finally, we list WFC3 calibration uncertainties. As noted in Table [H we find 0.01 magnitudes of tension with 
the STScI zeropoints. Our zeropoints also scatter by a few mmags depending on aperture radius, and show some 
tension between standard stars. Until these issues are resolved, we take a 0.01 uncertainty in the WFC3 bright 
zeropoints. Going from bright to faint zerop oints adds about 0.01 magnitudes of uncertainty for the WFC3 CRNL 
(|Riessl[201^ iRiess fc Petrol I 2 OIOI : lRiessll201ll) . and moves the effective zeropoints 0.04 magnit udes br ighter (lower). 
To be conservative, we also take half of our proposed update of the WFC3 bandpasses ('Section fA.3.2l as uncertainty, 
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giving 9 mmags in FllOW and 7 mmags in F160W. In total, we estimate that the WFC3 low-count-rate ST zeropoints 
are 28.428 dz 0.017 for FllOW and 28.176 dz 0.016 for F160W. Note that these zeropoints are tied to the CALSPEC 
system, which has uncertainties as well. 
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